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thresholds, and evaluation of the post— stability fluid— induced rotor lateral limit cycle 
self— excited vibration parameters (whirl and whip). 

Using the multi-mode modal approach [ 8 , 17], the linearised rotor system eigenvalue problem is 
solved, and approximate values of natural frequencies and thresholds of stability are given. It is 
shown that there exist fluid— generated natural frequencies of the system, and the corresponding 
sequence of modes. The nonlinear model allows for calculation of the after— threshold, self 
excited vibration parameters, including modes of vibration. 

2. MATHEMATICAL MODEL OF A TWO-MODE ROTOR WITH ONE SOURCE OF 
INSTABILITY 

Two physical models of rotors on which mathematical models are the same will be discussed in 
this section. In both cases the rotors carry one source of fluid-induced instability. 

Consider a balanced isotropic rotor of a fluid— handling machine supported by two rigid bearings. 
The rotor rotates concentrically in a seal located at the axial section "S" (Fig. la). Consider a 
b alan ced isotropic rotor rotating concentrically in one fluid— lubricated and one rigid bearing 
(Fig. lb). When using the multi— mode modal approach (mode summation), the mathematical 
model of the rotor, in both cases, is as follows: 

Mja j + Dji j + (Ki+K2)*i — Kjsj = 0 

Ms's 2 + D 2 S 2 + (K2+Ks)z2-K2*i-K 3 z3 = 0 

D(»3— jAflsj) + (K|j-fK3+K4)z3 + *3^*0 *3 1 ) — K 3 *2 = 0 

* i(t) = xi(t) + jyi(t), 1=1, 2, 3, I * 3 1 = V x 3 + y\ • i=V^T . * = d / dt 

where xj, yi are rotor horizontal and vertical displacements, Mj, M 2 , Kj, K 2 , K 3 , K 4 , Dj, D 2 are 
two-mode rotor modal masses, stiffnesses, and external damping coefficients, respectively (in case 
of the fluid-lubricated bearing, the stiffness K 4 either equals zero or represents an external 
supporting spring (Fig. lb)). D, K B , A are seal or bearing fluid radial damping, radial stiffness 
and circumferential average velocity ratio, respectively, Ip is fluid nonlinear stiffness function of 
radial displacement, 0 is rotative speed, t is time. 

For clarity of presentation, the unbalance forces, radial side— load forces, external cross damping, 
fluid inertia effect, and other nonlinear functions are omitted in the model. 
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2.1 


Ei genvalue Problem: Natural Fr e quencies. Thresholds of Stability 


r °' the Un “ ri " d E,> ' (i) (th “ »■ *>>' fo»owi»g 

D — -L _i_ V- i ~w 1 rrOr. 

( 2 ) 


where 


(«l« 2 -Kj)[jD(a>-Afl) + K B + Ks + K 4 ] - Kj«, = 0 
= Ki + Kj +1 + ju® j — MiW 2 , i=l ,2 


( 3 ) 

are rotor partial complex dynamic stiffnesses and u is the complex eigenvalue Eo (7) car, 
certainly be solved numerically for specific parameter values. R VS auAor’s ^ intentioT 
however, to show specific qualitative features of this equation solution. ’ 

The split of Eq. (2) into real and imaginary parts provides the following equations: 

DwC<^-An)(«,D2+« 2 D,) - («i« 2 -K^-u; 2 D 1 D 2 )(Kb+Kj+K 4) + K^/c, = 0 (4) 

D^Anj^Ks-K^DiD^ + w(/c 1 D 2 +« 2 Di)(Kb+K 3 +K 4 ) - K§o®, = 0 (5) 

stiffnesses ffl that “ real part, of the complex dynamic 

w Eo p 7r w-:^; 8 of thc systcm naturai frequende8 ^ ^ <* i 

provide, more accurate relationships: ^ ,2 ‘ Thc flrst a PP ro ™tion 

FLUID WHIRL FREQUENCY = a AQ Jl - 

— Di i^ 2+K 3 ~ M 2 A 2 fl 2 - f(Ki+K 2 -Mi A 2 fl 2 ) D 2 /Di] (Kn+K»+Ki)— K? ‘ 

D l (Kj +K 2 -M t A ^wi)(K 2 +K 3 -M 2 A W a/- k| - A 2 ii5D,D 2 ]J ( 6 ) 

FLUID WHIP FREQUENCIES = U> 2)3 « _ 

- (-llip, LK2+K ? -M 2 q, 2 ^(K, + K 2 -M , u; n 2 1)0,70 1 lfKp+K^-i-K ^l—y? . _ ^ 

2D(w n i — AW ) M (jlj) 4 > 1—1,2 (7) 

REVERSE MODE FREQUENCIES = ui 4l , a 

where ^,,^ 2 are solutions of the rotor partial characteristic equation /C,/C 2 - K] = 0 : 


Wni = 


^2 + Ki+K / [ Kj+ K 2 K 2 +K 3 ] 2 K l 

2M * 2M ’ 1 W [ 2M 1 - -2MTJ + wrk 


1/2 


, 1 = 1,2 ( 9 ) 


2^2^ w. “ nn °‘ b ' * PPli ' d f0t lh ' V,Ju ' «* '“« should he 

The fluid whirl frequency w, is almost a linear function of the rotative sneed TK. <1 m u- 
frequencies independent from the rntntiv, .peed. luhe'^X. “ P 


279 


in the plane (Q, a7i), i=l t 2,3, the straight lines degenerate into hyperbolas (Fig. 2), resulting in a 
smooth transitions from the whirl— to— whip and whip— to— whirl modes. 



Fig. 2 Natural frequencies versus rotative speed. 

Eq. ( 4 ) provides the approximate values for stability thresholds. For the whirl onset, U) in Eq. 
( 4 ) is substituted by Aflg T (^e zero— th approximation ( 6 ) with Di=D 2 =0). The resulting 

quadratic equation for establishes: 

K t -fK 2 K 2 , K 3 (K B +K 4 ) 

2M j 2M 2 2 M 2 (K 3 +KB+K 4 ) 


n ST i ~ X 




CKi+K 2 k 2 K 3 (Kb+K 4 ) 1 

3 . k$ 

[ 2M 1 2M 2 2M 2 (K s +Kg+K4)J 

+ MjM 2 


1/2 


i=l ,2 ( 10 ) 


There exist, therefore, two instability onsets leading to the fluid whirl instability. For the whip 
cessations, U) in Eq. (4) is substituted by W n j, 1=1,2 (the zero-th approximation (7) with Di = 
D 2 = 0), and f) is substituted by ftsTi+ 2 : 


n 


1 

ST i +2 “ X 


K? 

^ i+ DW n i[D 2 +D 1 (K 2 +K3-M2Wni 2 ) a /K!il 


i— 1,2 


( 11 ) 


Note that all stability thresholds are inversely proportional to the fluid circumferential velocity 
ratio A . The onsets ( 10 ) do not depend on damping, while the cessations ( 11 ) do: larger 

damping, lowers the cessation rotative speed, thus narrowing the range where the instability 
occurs. Practically, the cessations seldom occur at rotative speeds (11), because the rotor is 
already in the self— excited vibration condition, and the system parameters differ from the linear 
ones. (The first order approximation (11) does not show dependence on Kg, thus no 

dependence on if ) ; it is necessary to consider higher approximations). 

In the rotor/seal system case, the fluid radial stiffness Kg is usually smaller than the stiffness K 4 , 
thus the third term in Eq. (10), and the third term under the radical in Eq. (10), are 
approximately the stiffnesses of the rotor right portion, combined by two elements in sequence: 
KjK 4 /(Kj+K 4 ) (Fig- 1 & )- The expressions in large parentheses of Eq. (10) represent, therefore, 
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the natural frequencies of the entire shaft. The stability onset expressed as the 1/A fraction of 
the natural frequency (which for the most often used value A=0.5 equals 2) is well known in 
rotor/seal dynamics. At the rotative speed as high as double the first natural frequency, only 
whip type of instability may occur. That is why very often only whip, and no whirl vibrations, 
are associated to rotor/seal dynamics. As Eq. (10) indicates, however, the multiplier 1/A is not 
necessarily equal to 2 (which it very seldom is), but depends on the actual circumferential flow 
conditions in the seal. Swirl brakes or anti— swirl injections can significantly modify A [11]. 
Additionally, the fluid radial stiffness may modify the instability onset (10): the larger Kg, the 
higher the instability onset. (A larger fluid radial stiffness results, for instance, from a higher 
external pressure of the fluid.) The onset ( 10 ) also depends on the location of the seal along the 
rotor, thus on the values of stiffnesses K3 and K4* The highest values of the onsets ( 10 ) occur 
where the stiffness K3 is equal to Kg + K4, that is, when the seal is located approximately in the 
middle between the second modal mass and the right-side bearing. 

The function fj = K\K2 — K$, which provides the whip frequencies CJni, Eq. (9), and the 

functions f 2 — (/C1/C2— K\)/K\^ {$ = K§/(Kg+K3+K4), which provide onsets of instability, Eq. 
( 10 ), can be investigated graphically in the plane (u/^, f^), 1/= 1 , 2 , 3 . The result is presented in 
Figure 3 . The instability onsets can be found on intersections of the hyperbola f 2 and the 
straight line f 3 . Now it is easy to discuss the effect of the stiffnesses K 3i K4, Kg on the 
instability onsets. If the stiffness K4 is low (rotor/bearing model), then the line f 3 is higher, and 
the resulting stability onsets occur at lower rotative speeds, thus, very often in the whirl 
instability range. An increase of the stiffness K 3 (mass M 2 closer to the instability source) 
results, again, in a decrease of the instability onset. An increase of the fluid radial stiffness Kg 
provides a stability improvement by moving the instability onsets toward higher rotative speeds. 
In all cases, however, the onsets occur at the speeds fl gTi which are lower than Ufoi/A, i=l,2. 



Fig. 3 Graphical solution of Eqs. ( 4 ) and ( 5 ) showing whip frequencies ( 9 ) and onsets of 

instability (10). 


In this study the seal location is chosen outside the modal masses; a similar analysis can be 
performed for a case where the seal is located between the modal masses, that is being closer to 
the second mode nodal point. 
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2.2 Rotor Modes 

The complex modal functions that are calculated from the linear part of Eqs. (1) are as follows: 


$U/ — li $2 v — 1 + 


K rfD i j h>-M i (jfl 
K 2 


x _ K 3 (K|-t-K;-j-Di j U/— M i U 

™ ~ K 2 [Dj(t»A-AW)+KB + K 3 +K 4 ] 


(12) 


Using the natural frequencies (6) and (7), the approximate modal functions at stability 
thresholds are obtained. The modal functions of the first and second mode whirls are as follows: 


*»=!, ♦>!»{- H)‘ J(> + jjh , in - hi k b+ kU ' i=1 ' 2 


where 


? 2 [k 2 + 


1- Ml _ K 3 Mi(K R -fK 4 ) ' 
M 2 K 2 M 2 (K b +K 3 +K 4 ) . 


The modal functions of the first and second mode whips are as follows: 


(13) 


in.! = i, in,* « i + , 

D((d n i~ AOg T i+2 ) 

4>3W2» — foi+gK; e j K B +K 3 +K 4 j-^2 (14) 

V(Kb+K 3 +K 4 )^ + D*(Wbi-All SXi+3 )3 

where are given by Eq. (7) and flg Ti+2 b Y Eq* (11). 


During the whirl of the first mode, the entire shaft vibrates in phase. During the second mode 
whirl, the modal masses Mj, M 2 are 180* out of phase (when damping is neglected). During the 
first and second mode whips, two rotor masses vibrate in phase, and 180* out of phase, 
respectively, but the seal location leads the vibration of the second mass by an angle which 
depends on the system parameters (Eq. (14)). 

2.3 Self— Excited Vibrations 

After the onsets of instability, the rotor self— excited vibrations occur. Eqs. (1) have exact 
periodic solutions of the circular form: 


,1 = A,e j(Wt+ai ) , I2 = A2 e j(aA+a2) , , 3 = A,^ 


(15) 


with frequency (or frequencies) U) , and corresponding amplitudes A j , and phases ttj , i=l,2 
relative to the phase of the fluid instability source location. Eqs. (15) describe rotor self— excited 
vibrations, known as whirl and whip. They occur as steady— state limit cycles after the onsets of 
instability and some transient process. 


The frequencies, amplitudes, and phases of Eqs. (15) can be calculated from the set of nonlinear 
algebraic equations obtained by substituting (15) into Eqs. (1), and dividing all terms by e 1 ^. 



Note that the argument of the only nonlinear function Ip becomes now A3. The final algebraic 
equations are as follows: 


(o>— A42) ( /c 1 Kg— Kg~ 
LJ{ K1D2+K2D l) 


- k b - k 3 - k 4 


A ie j<*l = _Ml A) , Aae* ttg = A, (17) 

K 1 K 2 - Ki 2 - KJ 

where K\, i=l,2 are dynamic stiffnesses (3) and their real parts, respectively, and lp~^ denotes 

the inverse function ip . The frequency u) in Eq. (16) should be substituted by the frequencies 
obtained from the equation which differs from the characteristic equation (2) only by one 
nonlinear term. The approximate values for the self— excited vibration frequencies are, therefore, 
very close to the fluid whirl and whip frequencies. In Eqs. (6) and (7) the fluid radial stiffness 
Kb should read [Kb+V^A)]. Since the second terms in Eqs. (6) and (7) are proportional to 
rotor damping, and thus are relatively small, the above— discussed nonlinear function adjustment 
will, therefore, bring very little change to the approximate frequencies (8), (7). The reverse mode 
frequencies (8) do not satisfy Eq. (16). 

For the whip frequencies U) = (J n i <> r = k*n2i the denominators of Eqs. (17) are controlled only 

by damping, thus, independently of the amplitude A3 , the rotor with masses Mi, M2 is in 
‘'resonant" conditions corresponding to its two lateral modes. 

It is easy to show that for the first mode whirl and whip frequencies, the phase angles 0t\ and 
<*2 are negative, thus the mass locations are lagging the instability source at the bearing or Beal. 
For the second mode, the phase lag also occurs additionally to the relative 180" phase change of 
the modal mass locations. The instability source phase leading of the whirl or whip vibrations is 
well recognized in the field, and has become the basis of the identification method for the 
instability source axial location on machinery trains [18]. 

3. MULTI-MODE ROTOR WITH ONE SOURCE OF INSTABILITY 

An analysis similar to that presented in the previous chapter can be performed for the 
multi— mode rotor (Fig. 4). The rotor characteristic equation now has the following form: 

[Dj(uJ— Aft) + K B + K n -i + KnJJt — Kn-i = 0 (18) 

where 
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and 


*i - Ki + K 1+ , + jo;Di - MM i=l n _ 2 

are rotor partial complex dynamic stiffnesses. The stiffness k . t . . 

»• *• to. ^ th< Hght (Pi/rr ° r 

■"O' Wj ' p,ovi<1 ' th ' approximation to *° 



Fig- 4 Multi— mode rotor/bearing/seal model. 

S slis :iiXtcie?;s it 1 s *, Aa The further — *■ 

of Eq (18), are proportional to the external’ damping ^ cak ulated from the imaginary part 
Wni/A, they are small. ^ and except for the regions when f2 = 

Similarly to Eq. (4), the real part of Eq. (18) provides the stability thresholds. 

- * to performed for 

speeds lower than 1/A fraction of natural f . ' tbc insta bihty onsets occur at rotative 

parameters, namely, A, K n _,, K n , and Kjj. reqUCnCle8 ’ and can be controlled by the similar 


4.1 


experimental results 

Experimental R i p r 


disks of the mass 4.637 * 10 -3 lb.sVin^inbLd) and uS? i f 2 / diameter shaft with tw ° 

was supported inboard by a relatively rigid si d W * u k ! b ’ 8 /,n ( outboard )- The rotor 

oil-lubricated, cylindrical bearing with 1 5" diameter b “ ah, “?* ai | d outboard by a 360° 

was T10 oil with 5 psi inlet pressure. The IC ^ F&dlaI cle «ance. The lubricant 

flexible coupling. Four orthogonal springs supported thfrTtf b j e, ^ ctric motor through a 
maintain concentric position of the journal inside ti, k. ” b ° ard cnd of the shaft > in order to 
amount of unbalance: * m8lde at rcst - The rotor carried some 

other axial locations of th^taR (Fi^” h at ‘ he fluid ~| ubricated bearing and at four 

well as from the KeyphasoJ® (on ce otr JL vlbr fonal data from these transducers, as 
computerised data acquisition and processing system. ^ ** transducer * were collected by the 

4 -2 Experimental Results 
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SEC /oS I ? M0DE WHIRL: n = 9200 rpm, w, = 4220 cpm 
(TRANSITION TO SECOND MODE WHIP) 


■l 



Fig. 5 Rotor .ketch and rotor modal centerline, at the fir.t and .econd mode whirl, the fint 
mode whip, and at the tran.ition from the .econd mode whirl to the .econd mode whip 
at .pecific rotative .peed.. The orbit, are filtered to whirl or whip frequence. w if i=l,2 
Orbit Kale: 5 mil/div. Cue 1: disk, apart. 
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frequencies. The artificial Keyphasor marker positions were used to plot rotor centerline 
deflections for particular modes (Figs. 5 and 6). Figures 7 to 10 present the spectrum cascade 
plots during the rotor start— up. They exhibit synchronous (lx) unbalance— related vibrations 
and high rotor subsynchronous vibrations of the first and second mode whirl and of the whip 
type. Moving the disks toward the rotor midspan (case 2) resulted in a lowering the first natural 
frequency, and increasing the second one as compared with case 1. The first onset of instability 
occurred higher in case 2, but the second one occurred at a lower rotative speed. In both 
considered cases, rotor first mode whip vibration amplitudes were much higher than in the whirl 
modes. In both cases, the whip cessations occurred simultaneously with the second mode whirl 
onsets, but this is not a rule for other systems [9]. 



FIRST MODE WHIRL: fi = 2500 rpm, u)\ = 1200 cpm 

\ ■ ’l : * ’i ; ■ 


FIRST MODE WHIP: fl = 4000 rpm, u >2 = 1400 cpm 



SEC OND M ODE WHIRL: 0 = 8700 RPM, w, = 4200 cpm 


Fig. 6 Rotor sketch and rotor modal centerlines at the first and second mode whirl and at whip 
at specific rotative speeds. The orbits are filtered to whirl or whip frequencies UJ\, i=l,2. 
Orbit scale: 5 mil/div. Case 2: disks close to mid-span. 
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The modal parameter ratios as functions of the square root of mass ratio calculated from Eqs. 
(20) are presented in Figure 11. Since the modal parameter ratios must be positive, there exist 
limited ranges of mass ratios which satisfy this condition. It can be seen in Figure 11 that for 
case 1: 1.92 < fi < 10.8, and for case 2: 1.98 < /i < 19.9. 

Perturbation testing should provide more reliable identification data [17]. 



NOTATION 


Ai,ai 

Amplitudes and phases of rotor self- 
excited vibrations, respectively 

K\,K\ 

Rotor complex dynamic stiffnesses 
and their real parts, respectively 

d,k b 

Fluid radial damping and radial 
stiffness, respectively 

X 

Fluid circumferential average 
velocity ratio 

■i.lsil 

Rotor complex lateral displacements 
and their absolute values 

/*=Mi/M 2 Modal mass ratio 

Rotor complex modal functions 

Di,Mi 

Rotor modal damping coefficients 

1> 

Fluid nonlinear stiffness function 


and modal masses, respectively 

u 

Complex eigenvalue; also 

Ki 

Rotor modal stiffnesses 


self-excited vibration frequency 

i,*' 

j=v c i 

Integers 

Ui 

Approximate natural frequencies of 
the rotor Bystem 

t 

Time 

Wni 

Rotor natural frequencies 

*i.yi 

Rotor horizontal and vertical 

n 

Rotative speed 


displacements, respectively 

tt STi 

Thresholds of stability 
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During whip conditions the whip frequency filtered orbits are elliptical, reflecting rotor system 
anisotropic features. During whirl conditions the orbits are circular (except the first-mode whirl 
in case 2), and much smaller than at whips (Figs. 5, 6). There exist higher harmonic and 
combinational components in the spectra (Figs. 7 to 10), especially when amplitudes of whip 
vibrations are high (Figs. 7, 8) causing involvement of geometric nonlinearities. 

The first instability onset in case 2 (Figs. 9, 10) occurs at relatively high rotative speed, in the 
range of the transition to whip rather than at whirl, as in case 1 (Figs. 7, 8). The phases and 
amplitudes of these self-excited vibrations change considerably just after the onset: first the 

horizontal amplitudes are large, then, at higher rotative speed, in whip conditions, the vertical 
amplitudes become dominant (Fig. 9). 

4.3 Identification of Rotor Modal Parameters From Experime ntal Results 

In this section a method of partial (incomplete) identification of the rotor two-mode modal 
parameters will be outlined. This method does not pretend to be general; the identification 
results show only that the modal modelling is feasible, simple, and reflects major dynamic 
feature* of the rotor behavior. 

The data for identification is taken from the start-up runs (Figs. 7 to 10): rotor natural 

frequencies of the first two modes, Wm, W n 2. and the onsets of instability for the first and second 
mode whirls, v fl ST From Eqs. (9) and (10), the rotor parameters are calculated as 

follows: 

^ = [b(tdnl + - b) - = bfi - — 

( 20 ) 


= Wni + ItJnJ - b 


K; Knd-K4 _ K 3 

M 2 ’ Mj " M 2 


k 3 /m 2 


-1 


tt>nl+Wn2 — ^ a (Osx l+^ST 2 ) 


where 


k l<^n2 ~ 1 ^ST 2 

Wnl + Wn2 “ A*(fl|«j t + ^ST2) 


M = 


Mi 

M 2 


Since there are only four equations, (9) and (10), and seven unknown parameters, the 
stiffneas-to-mass ratios arc expressed as functions of the modal mass ratio . 


The data for the cases 1 and 2 are given in Table 1, 

TABLE 1. Experimental Data 



tdnl 

U n2 

^STl 

fi ST2 

H 

rpm 

— 

Case 1 

1677 

5200 

2000 

8800 

0.48 

Case 2 

1400 

6300 



0.48 
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